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Abstract 

Magnetic fingerprints on the spectra of an electron interacting with 
a negatively charged ion in a parabolic quantum dot (QD), and of 
two interacting electrons in such a dot, are investigated via a new 
pseudoperturbative methodical proposal. The effects of ion - electron 
and electron - electron interactions on the spectra are studied. The 
effect of the symmetry of such problem is emphasized. Compared with 
those obtained by Zhu et al.[6], via a series solution, the results are 
found in excellent accord. Higher excited - states are also reported. 
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1 Introduction 



Quantum dots (QDs) are essentially quasi - zero - dimensional little islands, 
populated by two - dimensional ( flatland) electrons laterally confined by an 
artificial potential [1,2]. The variation of their spectra with magnetic field is 
often called magnetic fingerprints of QDs. Recent intensive experimental and 
theoretical investigations have been carried out to explore various spectral 
properties of QDs. Yet such studies are in their infancy and still expanding 
rapidly [1-15]. 

The ion - electron and electron - electron interactions in such QDs are 
of great importance [6,8-11]. The Hamiltonians (section 2-1) of which are 
known to belong to non - exactly soluble Hamiltonians. One has therefore 
to resort to approximation methods to study such systems. There has been, 
in the last few years, an increasing interest in the study of two - electron 
QDs in a magnetic field [6,12-15]. Among the several eligible and refiable 
methods for solving the multi - electron problem in QDs exist: exact numer- 
ical diagonalization [2,16], numerical simulations based on quantum Monte 
Carlo techniques [17], Hartree - Fock calculations [16, 18-20], a series solu- 
tion based on asymptotically physical grounds of the wave functions in the 
regions < r and r < oo [6,21], etc. It is therefore interesting to carry out 
systematic studies of the electronic structures in QDs populated by one and 
two electrons with and without magnetic fields. 

Recently, we have introduced a pseudoperturbative ( artificial in nature) 
shifted - I { I is the angular momentum quantum number) expansion tech- 
nique ( PSLET) to solve Schrodinger equation for states with arbitrary num- 
ber of nodal zeros. It simply consists of using l/F as a pseudoperturbation 
parameter, where I — I — P and is a suitable shift to be determined [22- 
27]. Encouraged by its very satisfactory performance ( documented in [22-27] 
through comparisons with direct numerical integration, quasi perturbative. 
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variational, Hill determinant, etc, methods), we use PSLET recipe, in section 
2-2, and study the spectral properties of an electron in a parabolic QD with 
a negatively charged ion ( impurity), and two electrons in such a QD with 
and without the magnetic field. 

In section 3 we compare our results with the ones obtained by Zhu et 
al.[6], via a series solution, and report on other states that have not been 
considered before, to the best of our knowledge. The effect of ion - electron 
( section 3-1) and electron - electron ( section 3-2) interactions are also dis- 
cussed. In the absence of these interactions we show that whilst the magnetic 
field completely lifts the confinement degeneracy, as it increases from zero, 
for the one electron case, it partially hfts it for the two - electron case. More- 
over, at infinite magnetic fields the energy levels cluster around Landau ones, 
inducing in effect Landau degeneracy. The inclusion of these interactions, on 
the other hand, completely lifts the confinement and Landau degeneracies 
and change levels ordering. Consequently, energy levels crossing occur and 
spin - singlet spin - triplet oscillations are manifested for two - electron QDs 
in a magnetic field. Section 4 is devoted for concluding remarks. 

2 Hamiltonian models and PSLET recipe 
2.1 Hamiltonian models 

To a very good approximation, the harmonic oscillator describes the lateral 
confinement of electrons in some experimentally created QDs. The motion in 
the z - direction ( the growth direction) is always frozen out into the lowest 
subband [6,28]. The confinement in this direction is assumed to be stronger 
than that in the xy-plane, and the dots, in effect, can be treated as two - 
dimensional thin discs [6,21]. Then the Hamiltonian of an electron in such 
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a parabolic QD with a negatively charged ion center is very well simulated by 

H^0 = -V^ + \ll + (1) 

and that of two electrons in the same QD by 

Ho--Vl-Vl + \^j{pi + pl) + -^ (2) 

4 \Pl — p2\ 

( in effective atomic units) with 7^ related to the confinement region ( 
hence to the quantum size) of the electrons in the QDs. The energy and 
length are given in effective Rydberg R* and effective Bohr radius a*, respec- 
tively. 

When a magnetic field B is applied perpendicular to the xy - plane, 
through the symmetric gauge A — {—y,x,0)B/2, the Hamiltonians in (1) 
and (2), respectively, read 

Hi = Hio + ^7' + lL„ (3) 

and 

H = Ho + |7' (Pi + pI)+ 7^.1 + 7^.2. (4) 

Where 7 = fiuJc/2R* with the cyclotron frequency ujci and 7^21 and 7^22 are 
the induced Zeeman terms. Obviously, Hamiltonian (4) is separable and can 
be recast as 

H^Hn + Hr, (5) 
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where 



(6) 



and 



8 r 



(7) 



which represent the center - of - mass (CM) and the relative - motion 
(RM) Hamiltonians, respectively. Here, = 7^ + 7^, R = {pi + P2)/2, 
Vfl = Vi + V2, r = pi- p2, Vr = (Vi - V2)/2, and Lzr = -id/d(j) and 
Lzr — —id /dip are the Z- and z- components of the angular momentum oper- 
ators in the CM and RM systems, respectively. Effectively, the problem is re- 
duced into two: a quasi - particle of the CM coordinate in a rescaled external 
field, and a quasi - particle of the RM coordinate in a rescaled external field 
and a rescaled repulsive Coulomb field (emerging from the electron - electron 
interaction) . The cylindrical symmetry of the attendant problems invites the 
separability of the wave functions to obtain. Hence, the two - particle wave 
function in the plane polar coordinate systems is ^ K,M{R)^k,m{'i')e^^'^ e^^'^ ■ 
The Pauh principle demands that ^k,m{f)&^^'^ is antisymmetric with respect 
to inversion r — > — r. No restrictions on ^K,M{R)e^^'^ are imposed. One 
would therefore obtain spin - singlet (s = 0) and spin - triplet (s = 1) states 
for even and odd m, respectively, through the prescription s — (1 — (— 1)"*)/2. 
Eventually, Schrodinger equation for Hamiltonian (7) reads 



+ 



m 



1 



2„2 



Uk,ra{r) 



E{k, m) 



Uk,m{r), 



(8) 
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where Uk,m{f) = v^^fc.ml'"), and k and m are the radial and azimuthal 
quantum numbers in the RM system, respectively. 

It is convenient to define the electron - electron interaction energies as 

E^^{k,m) ^ E{k,m) - EoiKm), (9) 



where 

Eo{k,m) = {2 k+\m\ + l)T + m-f (10) 



are the energies of (8) without the Coulomb interaction. Then the total en- 
ergies of (4) are 

E(k,m;K,M) ^ E(k,m) + E(K,M), (11) 



where 

E{K,M) = {2 K +\M\ + 1)T + M-f (12) 



are the energies of (6), with K and M the radial and azimuthal quantum 
numbers in the CM system, respectively. 

Similarly, the ion - electron interaction energies are defined as 

Eie{k, m) = Ei{k, m) - Eio{k, m), (13) 



where 
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Eio{k,m) = {2 k+\m\ + l) T + m-f (14) 

are the energies of (3) without the ion - electron interaction and Ei{k, m) are 
the eigenenergies of in (3). 

2.2 PSLET recipe 

In the underlying radical radial time - independent Schrodinger equation, in 
h — m — 1 units, 



the isomorphism between orbital angular momentum / and dimensionality 
D invites interdimensional degeneracies to obtain. Which, in effect, adds a 
delicate nature to this equation and allows us to generate the ladder of excited 
states for any given k ( number of the nodal zeros in the wave function) and 
nonzero / from the / = result, with that A:, by the transcription D — > 
D + 21. For more details the reader may refer to rcf.s [26,29-31]. 

Our recipe starts with shifting the angular momentum quantum number 
I in (15) through I — I ~ [3 and use 1//^ as a pseudoperturbation expansion 
parameter. Hence, equation (15) reads 

(16) 

where Q is a constant that scales the potential V{q) at large - Id limit and 
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is set, for any specific choice of Id and k, equal to P at the end of the cal- 
culations. Here Id = I + {D — 3)/2, to incorporate the interdimensional 
degeneracies. Hence, I — > I = Id — P and D = 2 { with / = |m|) through 
out this paper. Next, we shift the origin of the coordinate system through 
X = l^^'^iq — qo)/qo, where Qo is currently an arbitrary point to be determined 
below. Expansions about this point, x — (i.e. q — Qo) would lead to 



1 d"^ a"^ ~ 



(17) 



with 



'fV{x{q)) = qil 



I 



1 ^V{q,) 



H Q 



+ l^'^B^x + Y.v^^\x)l-^'\ 



n=0 



(18) 



where 



(19) 



v^^\x)^-{2(5 + l)x + B^x\ 



(20) 



v^^\x) = S„+2 2:"+^ + (-l)"(2/3 + l) 



n+2 



(n + 1) 



X" 



+ ^(^_ti) _ 1^ . n>2, (21) 



^^^^ ^ ^-+i^^J;i!g' 



(22) 
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It is then convenient to expand Ek^i as 

oo 

Ek,i = E 4? (23) 

n=-2 

Equation (17), along with (18)-(22), is evidently the one - dimensional Schrodinger 
equation for a harmonic oscillator Q^x^/2, with = 2^2, and the remain- 
ing terms in Eq.(17) as infinite power series perturbations to the harmonic 
oscillator. One would then imply that 



(24) 



(-1) 



'k,l 



Where Qo is chosen to minimize Ef,^i , i. e. 

y = and — -4^ > 0. 
dqo dq^ 

Equation (26) in turn gives, with I — \/Q^ 

lD-P^^lqlV'{qo). 



(25) 



(26) 



(27) 



The shifting parameter j3 is determined by choosing /£^^ /'*=0. Hence 









= - 







3 + 



V'{qo) 



(28) 
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where primes of V{qo) denote derivatives with respect to Qo- Then equation 
(17) reduces to 



•1 oo 



n=0 



.n=l 



(29) 



Setting the wave functions with any number of nodes k as 



= Pk,i[x) exp{Uk,i{x)), 



(30) 



equation (29) readily transforms into the following Riccati equation: 



-1 oo 



n=0 



-E^o^r^^r" 



n=l 



(31) 



where the primes denote derivatives with respect to x. It is evident that this 
equation admits solution of the form 



n=0 



n=0 



(32) 



oo k—1 



F,,(x) = x^ + EE«S 

n=0 p=0 



(33) 



where 
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n+1 



m=0 



2m- 1 



0, 



(34) 



n+1 

E 

m=0 



2m 



(35) 



Substituting equations (32) - (35) into equation (29) implies 



- E [Uk^ + 4"^ r("+^)/2 ) 

2 n=0 ^ ^ 

{m)^x(n-m)7_n/2 , ^M^(n-m)7-(n+2)/2 



-1 oo n 



n=0 m=0 



n=0 



n=l 



CO 

^|^^(n)|--n/2^^(n)^-(n+l)/2^ 



.n=0 



(36) 



The solution of equation (36) follows from the uniqueness of power series rep- 
resentation. Therefore, for a given k we equate the coefficients of the same 
powers of I and x, respectively. One can then calculate the energy eigenvalues 
and eigenfunctions from the knowledge of Cm,n,k, Dm,n,k, and a^^j. in a hier- 
archical manner. Nevertheless, the procedure just described is suitable for 
a software package such as MAPLE to determine the energy eigenvalue and 
eigenfunction corrections up to any order of the pseudoperturbation series 
(23). 

Although the energy series, equation (23), could appear divergent, or, at 
best, asymptotic for small I, one can still calculate the eigenenergies to a very 
good accuracy by forming the sophisticated Fade approximation [22-24,26,32] 
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P^{1/1) = (Po + Pi// + • • • + PmII'')/{1 + qi/l + ■■■ + Qn/I^) (37) 



to the energy series (23). The energy series is calculated up to Ej^p/l^^ by 
E,, = + 45 + . . . + E^^f/l^' + 0{1/V\ (38) 

and with the Piq{1/1) Pade approximant it becomes 



E,,i[10,9]^FEi-'^ + P!,{l/l) 



(39) 



Our recipe is therefore well prescribed. 



3 Results and discussion 



It is obvious, to a scaling factor, that Hamiltonians (3) and (7) bear the same 
form of a hybrid of a repulsive Coulomb and oscillator potentials 



V{q) = b'q'/2 + 2/q. 



(40) 



Hence, equation (28) yields 



462 g3 _ 2 



(41) 



and, in turn, equation (27) reads 



y/b' qt - 2 qo- (42) 
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Once Qo is determined (often numerically) the coefficients Cm,n,k, Dm,n,k, 
and obtained in a sequential manner. Then the eigenvalues and 

eigenfunctions are calculated in one batch for each value of k, D, I, and b. 

In order to make remediable analysis of our results we have calculated 
the first twenty terms of our energy series. We have also computed the 
corresponding sequence of Pade approximants P2{1/1), P2{1/1),- ■ Pg{l/l), 
-^10(1/0 observed their effects on the leading energy term PE^f' . More- 
over, the twofold effect of the first term in the effective potential 

VefM = ^^^^^ + V{q) (43) 

should be in point. That is, for \m\ = it represents an attractive core that 
strengthens the comfinement, whereas for \m\ > 1 it represents a repulsive 
core and strengthens the Coulomb repulsion. 

3.1 Ion - electron interaction effect 

In table 1 we display the energies of an electron in parabolic QDs, including 
ion - electron interaction, in a magnetic field for 7^ = 0.2. They compare 
excellently with those reported by Zhu et al. in figures 1(a) and 1(b) of [6]. 
In addition we report PSLET results for 4d~, 4p~, 5d~, and 6f~ states. 

In figure 1 we plot the energies of an electron in such QDs ( with 7^ = 0.2) 
versus 7, excluding ion - electron interaction. Clearly, it shows that whilst 
the magnetic field completely lifts the well known confinement degeneracy, 
(2k + \m\ + l)7(i, as 7 increases from zero, it eventually introduces Landau 
degeneracy as 7 — > 00. And, at 7 = 00 only s - states are feasible and de- 
generate. That is, a state {k\m) would cluster at an s-state (A;, 0) through 
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the prescription 

2{k' - k) = -\m' \ - m . (44) 

In the range of small 7 ( namely, < 7 < 0.4), one clearly observes that 
there are minima for states with negative m, whereas for states with positive 
m and s-states (m = 0) increase monotonically. This should be attributed 
to the following: (i) for m — 0, the effective potential does not support a 
minimum and the Zeeman term has no effect on these states, hence s-states 
increase monotonically. (ii) For \m\ > 1, the effective potential supports 
minima and the Zeeman term 777,7 refines them for negative value of 777 or 
shifts them up for positive m. However, in the range of large 7 the magnetic 
field ( in ~ r^g^) dominates over the Zeeman term and the leading term in 
the effective potential, hence all energy levels increase monotonically in this 
range of 7. As a result of (i) and (ii) one observes the energy levels crossings 
between different states in the range of small 7 in figure 1. 

The inclusion of the ion-electron interaction ( figures 2-4) significantly 
changes the spectral properties mentioned above. At 7 = ( see figures 2-4), 
one notices that all energies are shifted up. For example, the Is-states is 
shifted (in Ry* units) by 0.6162, 2p- and 2p+ by 0.4666, Sd" and 3d+ by 
0.3829, 4f~ and 4f+ by 0.3305, etc. This shift decreases as |?77| increases. 
For higher levels ( larger |m|) g increases and the Coulomb repulsion ~ 1/g 
decreases, a characteristic which is reflected upon the energy shift - ups in 
figures 2-4. As 7 increases from zero, figures 2-4 clearly show that the ion- 
electron interaction completely lifts Landau degeneracy and results in change 
in levels ordering. Consequently, energy crossings occur and odd-even parity 
oscillations are manifested ( 2p+ crosses with 4d'~, 4d~ crosses with 4f+, 3d+ 
in figure 3). Figure 1(b) of Zhu et al.[6] also bear this out. 

To explain energy crossings, we consider the energy shift-ups as a result 
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of the ion-electron interaction. At 7 = 0.4 ( energy crossings occur in the 
range < 7 < 0.4) for example, the Is-state is shifted ( in Ry* units) by 
~ 1.03, 2p- and 2p+ by ~ 0.73, 3d- and 3d+ by ~ 0.59, 4f- and 4f+ by 
~ 0.50, and so on. Obviously the ion-electron interaction does not distinguish 
between positive and negative m values. Moreover, the shift-ups decrease as 
\m\ increase for a given k, or for a given \m\ as k increases ( because of 
the Coulomb repulsion characteristic mentioned above). This in fact explains 
why the energy ladder of the lower batch, say, of energy levels in figure 1 is 
turned up-side-down in figure 2. Similar scenarios could be also developed 
for figures 3 and 4. 

3.2 Electron - electron interaction effect 

Next, we calculate the e" - e" interaction energies Eee{nr, m) for two electrons 
in QDs and compare them, in tables 2 and 3, with those reported by Zhu 
et al. [6]. They are in almost exact accord. We also display PSLET results 
for 4s, 4p, 5d, 4f, and 5g states. Figure 5 shows that Eee{nr,m) increases 
with 7 and the levels ordering is -Eee(0, 0) > -Bee(l, 0) > -Eee(2, 0) > • • • > 
E^eiO, 1) > E^eih 1) > ^ee(2, 1) > • • • > E^^iO, 2) > E^^Ci-, 2) > • • •, etc. 

Before we proceed any further let us study the quantum size effect on the 
two-electron spectra. Table 4 shows the quantum levels of two electrons in 
QDs with different values of 7^(7^^''^)) hence showing the quantum size ef- 
fect, when 7 = 0. It documents changes in the levels ordering ( consequently, 
energy crossings and singlet-triplet spin oscillations occur) as 7^(7^ ^^^) de- 
creases ( increases) from 1(1). Hereby, it should be noted that Zhu's [6] result 
E{1, 0; 0, 1; 0) for 7^ = 0.2 is now corrected from 1,4402 to 1.4413 in accor- 
dance with his finding in table 1 of [6]. The levels orderings for — 0.05 
are changed from (h) (e) to (e) (h) and from (k) (n) to (n) (k) in accordance 
with our findings, based on the stability of, at least, the last five terms of 
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the Pade sequence ( a signal of convergence to the exact results) mentioned 
above. 

To study the e~-e~ interaction effect on the spectra, we plot the energy 
levels E{k,m; K, M; s) of two-electron QDs with 7^ — 0.2 and at different 
values of 7 are plotted in figure 6 ( excluding the e~ - e~ interaction) and 
figure 7 ( including e~ - e~ interaction). Figure 6 clearly shows that the 
magnetic field partially lifts the well known confinement degeneracy as 7 
increases from zero and induces Landau degeneracy as 7 — > 00. Again, one 
observes the clustering of the quantum levels around E{k, 0; K, 0; 0) states 
as 7 — > 00, following a similar trend as that of (44). 

The inclusion of the e~ - e~ interaction completely lifts the confinement 
and Landau degeneracies and changes the levels ordering. Consequently, 
energy crossings occur, at which singlet - triplet spin oscillations obtain, at 
specific magnetic fields. Figures 6 and 7 bear these out. One could follow a 
similar scenario as that in section 3-1 to explain the energy crossings, since 
the e~-e~ interaction is effectively a repulsive Coulomb term. 

In table 5 we display our results for two interacting electrons in QDs in a 
magnetic field at different values of 7, and for — 0.2. When plotted with 
the same scale, our results excellently agree with those of Zhu, figure 2(b) in 
[6]. However, it should be mentioned that the quantum levels classified in 
figure 2(b) of Zhu [6] as a, b, c, ■ ■ ■ do not bear positive values of m and M 
but rather negative ones, classified in table 5 of this text. 

4 Concluding remarks 

In this paper we have used our recently developed PSLET theory [22-27] to 
study the magnetic fingerprints on the spectra of an electron in parabolic QDs 
with negatively charged ion, and on two interacting electrons in such QDs. 
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We have also emphasized the effect of the symmetry of the problem ( marked 
in the leading term of the effective potential in (43)) on the confinement and 
Coulomb repulsive terms. The comparison of PSLET results with those of 
Zhu et al.[6] is readily very satisfactory. 

Although we have started with the central force Schrodinger equation 
and augmented the orbital angular momentum by I — > Id — I + {D — 3)/2, 
to incorporate interdimensional degeneracies, we have only considered the 
D — 2 with Z = |m| case for the attendant problems. 

A general observation concerning the method used by Zhu et al.[6] is in 
order. We have already mentioned that the series solution method used by 
Zhu et al.[6] is based upon the asymptotically physical grounds of the wave 
functions in the regions < r and r < oo ( i.e. r — > and r — > oo, 
respectively). Effectively and obviously, the authors used the asymptotic 
behaviours of the wave functions at weak and strong magnetic field lim- 
its implicitly. In the weak field limit the wave function is dominated by a 
Coulombic character ( hence Coulombic hke basis are used in equation (12) 
of [6]) and in the strong field limit its dominated by a harmonic oscillator 
character ( hence harmonic oscillator like basis are used in equation (14) of 
[6]). In fact, this is the only explanation, we could think of, as to why our 
results do not exactly agree with those of Zhu in table 2, for F = 0.4. 

The conceptual soundness of our PSLET is obvious. It avoids troublesome 
questions pertaining to the nature of small parameter expansions ( weak or 
strong field limits), and trend of convergence to the exact results. On the 
computational and practical methodical sides, it offers ( beyond its promise 
as being quite handy) a useful pseudoperturbation prescription where the 
zeroth order approximation P E^i'^^ inherits a substantial amount ( more 
than 90% for the above problems) of the total energy. 

Finally, PSLET theory could be applied to two electrons in an external 
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oscillator potential in D = 3 space [33], QD lattices [34], 3 - electron QDs 
[35], 2 - dimensional hydrogenic donor states in a magnetic field [36-39], etc. 
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Figures captions 

Fig.l: Ei{k,m) versus 7 for the ion - electron in QDs with 7^ = 0.2, 
excluding the ion - electron interaction. 

Fig. 2: Same as figure 1 including the ion - electron interaction for Is, 
2p~, 3d~, and 4f~ states. 

Fig.3: Same as figure 2 for 2s, Sd"*", 4f+, 3p~, 2p+, and 4d~ states. 

Fig.4: Same as figure 2 for 3s, Bp"*", 4p~, 5d~, and 6f~ states. 

Fig.5: Electron - electron interaction energies Ef.e{k,m) versus F for the 
states reported in tables 2 and 3. 

Fig.6: E{k, m; K, M) versus 7 for two electrons in QDs with 7^ = 0.2 
excluding the e~ - e~ interaction. 

Fig. 7: Same as figure 6 including the e~ - e~ interaction. 
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Tabic 1: Ei[k,m) with the ion - electron interaction for 7^ = 0.2 and 
different values of 7 



/ 





0.1 


0.2 


0.3 


0.4 


EAO 0) (Is) 


0.8162 


8860 


1.0537 


1.2617 


1.4816 


Ei{0, -1) {2p-) 


0.8666 


0.8445 


0.9339 


1.0724 


1.2282 


Ei{0, -2) {3d-) 


0.9829 


0.8775 


0.9100 


1.0070 


1.1307 


E,(0,-3) (4/-) 


1.1305 


0.9449 


0.9278 


0.9922 


1.0922 


E,iO, 1) (2p+) 


0.8666 


1.0445 


1.3339 


1.6724 


2.0288 


Ei{0,2) {3d+) 


0.9829 


1.2775 


1.7100 


2.2070 


2.7307 


M0,3) (4/+) 


1.1305 


1.5449 


2.1278 


2.7922 


3.4922 


^•.(1.0) (2s) 


1.1741 


1.2864 


1.5612 


1.9097 


2.2867 


E,(l,-1) (3p-) 


1.2336 


1.2557 


1.4563 


1.7415 


2.0623 


E,{l,-2) {4d-) 


1.3599 


1.2999 


1.4466 


1.6940 


1.9860 







0.05 


0.1 


0.15 


0.2 


Ei{l,l) (3p+) 


1.2336 


1.3156 


1.4557 


1.6405 


1.8563 


£;,(2,2) (3s) 


1.5406 


1.5815 


1.6967 


1.8692 


2.0814 


E,{2,-1) (4p-) 


1.6077 


1.6014 


1.6748 


1.8101 


1.9888 


E,(2,-2) (5rf-) 


1.7415 


1.6900 


1.7273 


1.8339 


1.9893 


E,(2,-3) (6/-) 


1.9002 


1.8042 


1.8069 


1.8869 


2.0218 
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Table 2: Electron - electron interaction energies Eee{k, m) at different values 
of r. The values in brackets are reported by Zhu et al.[6]. 



r 




0.05 


0.1 


0.2 


0.4 


£^ee(0, 0) 


(Is) 


0.1963 


0.3081 


0.4816 


0.7479 






(0.1963) 


(0.3081) 


(0.4816) 


(0.7494) 


Eee{l. 0) 


(2s) 


0.1853 


0.2871 


0.4413 


0.6716 






(0.1853) 


(0.2871) 


(0.4413) 


(0.6723) 


EeeM 


(3s) 


0.1763 


0.2703 


0.4106 


0.6164 






(0.1763) 


(0.2703) 


(0.4106) 


(0.6000) 


^ee(3,0) 


(4s) 


0.1686 


0.2565 


0.3863 


0.5742 


^ee(0, 1) 


(2p) 


0.1562 


0.2333 


0.3451 


0.5057 






(0.1562) 


(0.2333) 


(0.3451) 


(0.5066) 


^ee(l, 1) 


(3p) 


0.1468 


0.2168 


0.3170 


0.4597 






(0.1468) 


(0.2168) 


(0.3170) 


(0.4637) 


See (2, 1) 


(4p) 


0.1392 


0.2039 


0.2959 


0.4265 


£^ee(0,2) 


(3d) 


0.1311 


0.1915 


0.2776 


0.4000 






(0.1311) 


(0.1915) 


(0.2776) 


(0.3998) 


See(l,2) 


(4d) 


0.1240 


0.1799 


0.2594 


0.3720 






(0.1240) 


(0.1799) 


(0.2594) 


(0.3784) 


-Eee(2, 2) 


(5d) 


0.1182 


0.1707 


0.2452 


0.3505 


See(0,3) 


(4f) 


0.1144 


0.1652 


0.2375 


0.3399 




(5g) 


0.1025 


0.1472 


0.2105 


0.3002 
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Table 3: Same as table 2 for different values of F. 



i 




1 

i 


9 ^ 
Z.O 


A 


r. 





{Is) 


1 Q 1 Ofi 

i.oiyo 




9 nOQ/l 


Q Qono 

o.oyyz 








('9 980*^ 


('9 QQ'^n'l 






[ZS J 


11/1 7Q 


1 01 


9 /I 791 


9 7SR9 






^ 1 . 1^ ( O J 




l'9 ZL791 






[OS ) 


1 DQQQ 


i.oyoo 


9 17/1/1 
Z. i ( 44 


o /I /I qcr 
Z.44O0 






(^i.Uooo ) 


(^i.oyo ( ) 


[Z. I 1 44 ) 


('9 /1/1 
(^Z.4400 ) 




(As) 


0.9511 


1.5469 


1.9754 


2.2166 


EeeiO, 1) 


(2p) 


0.8279 


1.3404 


1.1707 


1.9195 






(0.8279) 


(1.3404) 


(1.7107) 


(1.9195) 


£^ee(l, 1) 


(3p) 


0.7439 


1.1941 


1.5188 


1.7017 






(0.7438) 


(1.1941) 


(1.5188) 


(1.7017) 


Eeei^, 1) 


(4p) 


0.6860 


1.0964 


1.3922 


1.5589 


^ee(0,2) 


(3d) 


0.6436 


1.0294 


1.3077 


1.4645 






(0.6436) 


(1.0294) 


(1.3078) 


(1.4645) 


£^ee(l,2) 


(4d) 


0.5957 


0.9496 


1.2047 


1.3485 






(0.5957) 


(0.9496) 


(1.2047) 


(1.3485) 


^ee(2,2) 


(5d) 


0.5598 


0.8907 


1.1292 


1.2636 


^ee(0, 3) 


(4f) 


0.5432 


0.8650 


1.0969 


1.2277 


^ee(0,4) 


(5g) 


0.4782 


0.7598 


0.9628 


1.0772 
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Table 4: E{k, m; K, M; s) quantum levels of two electrons in QDs at different 
values of 7^(7^^^^) for 7 = 0, including electron - electron interaction. Zhu's 
results [6] are obtained by replacing the last j digits of our results with the j 
digits in parentheses. 



7<i(7d ^^^) 


1(1) 


0.4(1.5811) 


0.2(2.2361) 


0.05(4.4721) 


a:(0,0;0,0;0) 


a) 3.3196 


a) 1.5479 


a) 


0.8816 


a) 0.2963(2) 


b:(0,l;0,0;l) 


b) 3.8279(8) 


b) 1.7057 


b) 


0.9451(0) 


b) 0.3062 


c:(0,0;0,l;0) 


c) 4.3196 


c) 1.9479 


d) 


1.0776 


d) 0.3311(0) 


d:(0,2;0,0;0) 


d) 4.6436 


d) 2.0000 


c) 


1.0816 


c) 0.3463(2) 


e:(0,l;0,l;l) 


e) 4.8279(8) 


e) 2.1057 


e) 


1.1451(0) 


e) 0.3562 


f:(l,0;0,0;0) 


f) 5.1473(2) 


f) 2.2716 


h) 


1.2375(156) 


h) 0.3644(476) 


g:(0,0;l,0;0) 


g) 5.3196 


h) 2.3399 


f) 


1.2413(02) 


i) 0.3811(0) 


h:(0,3;0,0;l) 


h) 5.5432(174) 


g) 2.3479 


i) 


1.2776 


f) 0.3853(4) 


i:(0,2;0,l;0) 


i) 5.6436 


i) 2.4000 


g) 


1.2816 


g) 0.3963(2) 


j:(l,l;0,0;l) 


j) 5.7439(8) 


j) 2.4597 


j) 


1.3170 


j) 0.3968 


k:(0,l;l,0;l) 


k) 5.8279(8) 


k) 2.5057 


k) 


1.3451(0) 


n) 0.4025(66) 


1:(1,0;0,1;0) 


1) 6.1473(2) 


1) 2.6716 


n) 


1.4105(053) 


k) 0.4062 


m:(0,0;l,l;0) 


m) 6.3196 


n) 2.7002 


1) 


1.4413 


0) 0.4240 


n:(0,4;0,0;0) 


n) 6.4782(93) 


m) 2.7479 


0) 


1.4594 


p) 0.4311(0) 


o:(l,2;0,0;0) 


0) 6.5957(6) 


0) 2.7720 


P) 


1.4776 


1) 0.4353(4) 


p:(0,2;l,0;0) 


p) 6.6436 


p) 2.8000 




1.4816 


m) 0.4463(2) 
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Tabic 5: E^kjin; K, M; s) quantum levels of two electrons in QDs with 
different values of 7 and for 7^ = 0.2, including electron - electron interaction. 



7 





0.05 


0.1 




0.2 


0.3 


0.4 


A:(0,0;0,0;0) 


0.8816 


0.9033 


0.9644 


1 


1664 


1.4217 


1.6967 


B:(0,-1;0,0;1) 


0.9451 


0.9194 


0.9380 


1 


0667 


1.2599 


1.4790 


C:(0,0;0,-1;0) 


1.0816 


1.0595 


1.0880 


1 


2493 


1.4822 


1.7439 


D:(0,-2;0,0;0) 


1.0776 


1.0067 


0.9889 


1 


0648 


1.2211 


1.4129 


E:(0,-1;0,-1;1) 


1.1451 


1.0756 


1.0616 


1 


1496 


1.3201 


1.5262 


F:(0,-3;0,0;0) 


1.2375 


1.1220 


1.0697 


1 


0984 


1.2249 


1.3960 


G:(1,0;0,0;0) 


1.2413 


1.2742 


1.3670 


1 


6766 


2.0734 


2.5066 


H:(1,-1;0,0;1) 


1.3170 


1.3030 


1.3548 


1 


5963 


1.9378 


2.3238 


I:(l,0;Orl;0) 


1.4413 


1.4304 


1.4906 


1 


7594 


2.1339 


2.5538 


J:(lr2;0,0;0) 


1.4594 


1.4004 


1.4166 


1 


6078 


1.9158 


2.2773 


K:(1,-1;0,-1;1) 


1.5170 


1.4592 


1.4784 


1 


6792 


1.9984 


2.3710 
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Fig.3 




Fig.7 
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